Improved models of upper-level 
wind for several astronomical 
observatories 



Lewis C. Roberts, Jr 1 * and L. William Bradford 2 

1 Jet Propulsion Laboratory, California Institute of Technology 
4800 Oak Grove Dr, Pasadena CA 91009, USA 

^Pacific Defense Solutions, 
1300 N. Holopono St. Suite 116, Kihei, HI 96753, USA 



*lewis.c. roberts@jpl.nasa.gov 
Accepted to Optics Express 

Abstract: An understanding of wind speed and direction 
as a function of height are critical to the proper modeling of 
atmospheric turbulence. We have used radiosonde data from 
launch sites near significant astronomical observatories and 
created mean profiles of wind speed and direction and have also 
computed Richardson number profiles. Using data from the last 
30 years, we extend the 1977 Greenwood wind profile to in- 
clude parameters that show seasonal variations and differences 
in location. The added information from our models is useful 
for the design of adaptive optics systems and other imaging 
systems. Our analysis of the Richardson number suggests that 
persistent turbulent layers may be inferred when low values 
are present in our long term averaged data. Knowledge of the 
presence of these layers may help with planning for adaptive 
optics and laser communications. 
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1. Introduction 

Images collected with ground based telescopes show temporal variation in im- 
age quality, primarily due to fluctuations in the index of refraction of the air 
along the light path. These index of refraction fluctuations are caused by differ- 
ences in air temperature arising from turbulent mixing. These turbulent motions 
happen on smaller scales than the gross motions of the winds. Still it is known 
that the wind has a part in the strength of these fluctuations. Hufnagel [Q]| de- 
rived a model for the vertical profile of strength of turbulence C„(z) which de- 
pended on a wind speed profile centered at a particular altitude. Furthermore, it 
has been observed that the index of refraction fluctuations are enhanced where 
there are shears in winds Q, so that direction as well as speed can be a factor, 
although this is difficult to quantify 0. 

Aside from the effects of wind on turbulence strength, the movement of tur- 
bulent regions past an observer induces temporal variations in seeing. If wind 
speeds are sufficiently high that turbulence does not evolve significantly over 
a measurement period, we apply Taylor's hypothesis and treat the temperature 
fluctuations as quantities that are advected by the wind flU. 

The following wind model is often used in the modeling of adaptive optics 
and other imaging systems (e.g. 0171 [SI), 



where v(z) is the wind speed in m/s, z is the height in m, and is the angle from 
the zenith. Bufton [5 J is often cited as the source for this model (e.g. (Hill El), 
but that paper does not present an explicit wind model. 

The basic form of the model was presented by Greenwood Readers often 
over look the statement in the body of the original paper that the height z = 
corresponds to a mean sea level altitude of 3048 m. This is the height of the 
Maui Space Surveillance System on Haleakala, where the analysis was being 
used. So the height of the tropopause term, 9400 m, should really be 12448 m 
when the model is used for sites other than Haleakala. This caveat has been 
overlooked in a number of references, leading to misleading results. 

The model from Greenwood [0 was derived by averaging radiosonde data 
collected from balloons launched from Lihue on the island of Kauai from 1950- 
1970 and from Hilo on the island of Hawaii from 1950-1974. This was more 
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fully detailed in a technical report, which did not explicitly generalize the equa- 
tion, but the authors clearly understood what the various terms signified lITOll . 
The wind model is generalized as, 



where v is the wind speed, z is the altitude, vq is the wind speed at the ground 
or low altitude, vj is the wind speed at tropopause, hj is the altitude of the 
tropopause, L T is the thickness of the tropopause layer and £ is the zenith angle 
of the observation ifTTTl . 

There are several uses for a wind profile. It can be used in conjunction with 
a turbulence model to predict the performance of instrumentation. Most com- 
monly this is applied to the servo bandwidth of closed loop systems [9|. The 
minimum bandwidth of the system (to limit servo lag) is related to the Green- 
wood frequency. Because the Greenwood frequency depends on a product of 
and a power of the wind speed v 5 / 3 , it is important to get the wind profile 
and turbulence profiles as a function of altitude. A wind profile can also be used 
to simulate turbulence in conjunction with phase screens. 

Currently, there is an increased emphasis on modeling effects of turbulence 
on telescopes with diameters of 20-50 meters. Significantly more spatial varia- 
tion of wavefront optical path difference is captured with such large telescopes 
in a short exposure than on a smaller telescope. Some of that variation is due 
to the presence of multiple layers of turbulence, moving at different speeds and 
in different directions. The patterns may not simply translate across the tele- 
scope aperture as implied by the Taylor hypothesis, but are more likely to "boil" 
with a complicated wind profile. This can have an impact on the design and 
performance of an adaptive optics system. 

In order to better model atmospheric turbulence, we used archival radiosonde 
data to compute the wind profiles for a number of astronomical observatories. 
We also show that the wind profile has a strong seasonal variation. We also 
extend the wind profiles, by including wind direction. Unfortunately, the wind 
directions do not always lend themselves to an analytic expression, so we only 
show them in graphs. 

Using the same radiosonde data, we examined the variation of the gradient 
Richardson number lfT2lO . which is an indicator of the stability of parcels of 
air. Since it depends upon strength of wind shear, as opposed to wind speed, 
it also serves to inform the modeler that wind speed alone is not sufficient. As 
will be seen in our plots, the wind speed can be modeled with a few simple 
terms, but the form of the Richardson number plots indicates the true nature 
of atmospheric turbulence: layered, with sharp boundaries (although these are 
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muted by the resolution of our data). 



2. Data Analysis 

Worldwide there are over 900 sites that launch radiosondes on a routine ba- 
sis. Radiosondes are small instrument packages mounted on weather balloons. 
They record atmospheric pressure, temperature, potential temperature, mixing 
ratio, dew point, and relative humidity as a function of height. Data are taken 
at about 6 second intervals, but databases may only have the data recorded at 
60 second intervals. Readings at certain specific standard barometric pressures 
are made regardless of when they occur. Since the balloons rise at about 5 m/s, 
the best vertical resolution would be approximately 30 meters, while the usual 
recorded resolution would be 300 meters and frequently the recorded data are 
at irregular intervals. Wind speed and wind direction are obtained from either 
radar observation or a navigation system such as LOng RAnge Navigation (LO- 
RAN) and newer systems may use the Global Positioning System. Almost all 
of these sites launch balloons at UT and 12 UT. The information is used as 
entries to weather prediction and simulation programs. 

We downloaded atmospheric sounding data for 1973 January till 2006 
September for a number of stations from the University of Wyoming's Weather 



Web (http://weather.uwyo.edu/upperair/sounding.html). The stations were cho 



sen based upon their proximity to an astronomical observatory. Not all the sites 
had data for the entire time period; the dates for the available data are also listed 
in the table. In addition there are gaps of varying length in each data set, which 
are most likely due to equipment problems. These data drop outs are not a sig- 
nificant problem due to the large number of data points that are available. The 
names of stations and their details along with the nearby observatories are listed 
in Table [Q The name of the station is given, along with its latitude, longitude 
and altitude, the dates when data were recorded and the names of nearby obser- 
vatories. In parenthesis are the distance in km between the observatory and the 
radiosonde launch site and the altitude of the observatory in km. 

For each sounding the pressure, temperature, dew point, relative humidity, 
wind speed and direction are measured as a function of height. The heights 
are irregularly measured and not repeated, so the data were interpolated onto a 
fixed height grid with an increment of 400 m and the mean was computed for 
various time periods. The height grid is measured from mean sea level, rather 
than the local surface. The different soundings have data measured to varying 
maximum altitudes, depending on when the balloons burst. We did not analyze 
any data from altitudes higher than 30 km, since there were too few data points 
for a meaningful analysis. 

Fig. [T] shows the yearly variations that occur for January and July in Hilo, 
Hawaii. Plotted over this are the mean wind speeds for the respective month 



Table 1. Details of the radiosonde launch sites used 
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for the entire 33 year period of data. This shows the variation in monthly mean 
wind speeds. Some of this difference is caused by large scale weather patterns 
and some is the natural chaotic behavior of weather. The overall shape of the 
wind profiles is highly consistent over the 33 years of data. 




Fig. 1 . The solid line in the figure on the left is the model fit to the mean 
monthly wind speed for Hilo, Hawaii in the month of January. The dot- 
ted lines are the mean wind speed for Hilo, Hawaii in each January from 
1973-2006. The figure on the right is the same, but for July. These figures 
illustrate the variation in the monthly mean wind speed. 



In addition to the sites listed in Table Q3 we also downloaded the wind pro- 
files for Lihue HI. As mentioned above this was used in Greenwood's original 
model. We computed the correlation coefficient between the Hilo and Lihue 



data sets for each month of the mean wind profiles. The correlation was highest 
in February (0.999) and lowest in August (0.981). The average correlation was 
0.992 with a standard deviation of 0.006. This correlation was high enough that 
we felt it was valid to only use the Hilo data. This does not say that Hilo and 
Lihue constantly experience the exact same wind conditions, but that on aver- 
age the wind profiles are very similar. This also shows that it is valid to use the 
same wind model for the entire state of Hawaii including the observatories on 
Haleakala and Mauna Kea. 

Surface winds are highly dependent on local geography, such as being in the 
lee of a mountain, and are also time dependent due to on-shore and off-shore 
breezes for the island sites and katabatic flows for sites near mountains. They 
are also dependent on the solar heating of the ground. This is less of a problem 
for astronomical uses of the wind profiles since mountain top observatories are 
usually situated above these surface effects. In addition, the wind profiles do 
not capture the effect that the local orography will have on the wind profile. 

We created the mean wind speed and direction for each of the 12 months 
for each of the sites listed in Table [Q We used a non-linear least squares fit 
algorithm to fit the data to 



Tables list the monthly coefficients and their associated error bars for all 
sites listed in Table CD The error bars are computed by the non-linear least 
squares fit. The errors bars increase as the winds vary from the pure Gaussian 
model. We fit the Gaussian to the wind speed data from the ground to a height 
of 16 km, beyond this the stratospheric winds start to affect the fit. The effect is 
marginal and often below the error bars, except for the summer data sets. 

The Gaussian wind model of |[9l (Eq. [U only fit the troposphere winds and 
ignored the stratospheric winds. It is possible to fit both of those winds if a 
Gaussian is summed with a second order polynomial. This has limited utility 
to astronomical imaging, as the turbulence above 20km is relatively weak, but 
high altitude airships are being designed to operate between 20 and 30 km [fT3l 
and require a knowledge of the wind speeds and directions in this region. 

The results are shown in Figs. [2]-[8] respectively. Each figure is composed of 
12 graphs corresponding to the 12 months of the year. The first letter of the 
corresponding month is in the upper right corner. Each graph has the the wind 
speed plotted with a solid line and the fitted wind speed computed from the 
coefficients in Tables [2J48] is plotted as a dashed line. The wind direction is 
plotted as a solid red line. The wind speed is marked along the bottom of the 
graph, while the wind direction is marked on the top of the graph. The wind 
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Table 2. Hawaii Wind Speed Model Coefficients 



Month 


AO 


A 1 

Al 


A <~\ 

A2 


A3 


Jan. 


4.1 ± 0.4 


25.8 ± 0.5 


12007 ± 69 


A f\ A 1 

4047 ± 1 15 


Feb. 


4.2 ± 0.4 


27.6 ± 0.4 


H952 ± 61 


3957 ± 99 


Mar. 


5.0 ±0.3 


29.6 ± 0.4 


12335 ± 55 


3405 ± 79 


Apr. 


5.0 ±0.3 


29.4 ± 0.4 


12776 ± 51 


3001 ± 68 


May. 


4.8 ±0.4 


22.9 ±0.5 


12909 ± 88 


3161 ± 117 


Jun. 


4.8 ±0.3 


17.0±0.5 


12637 ± 95 


2877 ± 127 


Jul. 


5.0 ±0.3 


14.0±0.5 


12406 ± 102 


2605 ± 133 


Aug. 


4.6 ±0.3 


12.8 ±0.5 


12527 ± 99 


2575 ± 128 


Sep. 


4.2 ±0.3 


13.6 ±0.4 


12263 ± 93 


2868 ± 126 


Oct. 


4.5 ±0.3 


17.2 ±0.4 


12222 ± 75 


2845 ± 101 


Nov. 


5.1 ±0.3 


17.4 ±0.5 


12313 ±92 


3214 ± 129 


Dec. 


5.0 ±0.4 


21. 1 ±0.5 


12199 ±84 


3617 ± 126 



Table 3. Oakland Wind Speed Model Coefficients 



Month 


AO 




Al 






A2 


A3 




Jan. 


-1.4 ± 1.8 


30, 


,5 


± 


1 


.7 


10730 ± 81 


8329 


± 


318 


Feb. 


-0.7 ± 1.6 


30, 


,0 


± 


1 


.5 


10695 ± 77 


8047 


± 


283 


Mar. 


-0.9 ± 1.4 


29, 


,5 


± 


1 


.3 


10548 ± 66 


7878 


± 


245 


Apr. 


0.6 ± 1.0 


28 


,3 


± 


1 


.0 


10472 ± 59 


7195 


± 


188 


May. 


2.8 ±0.7 


22, 


,8 


± 





.7 


10649 ± 70 


6367 


± 


170 


Jun. 


3.1 ±0.7 


19 


7 


± 





.7 


11 164 ±98 


6430 


± 


215 


Jul. 


4.0 ± 0.6 


16 


,5 


± 





.6 


11581 ± 122 


5540 


± 


210 


Aug. 


4.4 ± 0.5 


17 


,3 


± 





.6 


11602 ± 113 


5187 


± 


184 


Sep. 


3.2 ±0.6 


18 


,8 


± 





.6 


11400 ±95 


6240 


± 


192 


Oct. 


1.3 ±0.8 


22, 


,3 


± 





.8 


11073 ±79 


7297 


± 


217 


Nov. 


-2.0 ± 1.6 


32, 


,8 


± 


1 


.5 


10878 ± 72 


8251 


± 


263 


Dec. 


-1.8 ±2.0 


31 


,1 


± 


1 


.8 


10611 ± 81 


8358 


± 


333 



Table 4. San Diego Wind Speed Model Coefficients 



Month 


AO 


A 1 

Al 


A ^\ 

A2 


A3 


Jan. 


-2.9 ± 1.4 


33.9 ± 1.3 


11470 ± 76 


6104 ± 249 


Feb. 


1.4 ± 1.1 


32.4 ± 1.0 


11400 ± 81 


5081 ± 199 


Mar. 


-1.3 ± 1.2 


32.7 ± 1.1 


11701 ± 82 


5793 ± 226 


Apr. 


-2.4 ± 1.4 


32.7 ± 1.3 


11656 ±88 


6062 ± 266 


May. 


0.2 ± 1.2 


25.6 ± 1.1 


11577 ± 113 


5471 ±295 


Jun. 


3.6 ±0.9 


17.6 ±0.9 


12081 ±201 


4546 ± 364 


Jul. 


3.7 ±0.6 


10.6 ±0.6 


11483 ± 180 


3909 ±317 


Aug. 


4.4 ± 0.4 


10.5 ±0.5 


11821 ± 161 


3318 ±239 


Sep. 


4.4 ± 0.6 


18.8 ±0.6 


11897 ± 121 


3745 ± 193 


Oct. 


2.2 ±0.8 


22.7 ± 0.7 


11805 ± 111 


4843 ± 230 


Nov. 


-1.6 ± 1.3 


30.6 ± 1.2 


11511 ± 88 


5776 ± 256 


Dec. 


-1.3 ± 1.2 


32.2 ± 1.1 


11354 ±76 


5639 ± 223 



Table 5. Tucson Wind Speed Model Coefficients 



Month 


AO 


Al 


A2 


A3 




Jan. 


-1.1 ±0.7 


34.1 ±0.6 


11467 ± 34 


5549 


± 


112 


Feb. 


1.6 ±0.6 


34.1 ±0.6 


11620 ±40 


4993 


± 


104 


Mar. 


1.4 ±0.7 


34.0 ± 0.6 


11779 ±46 


5163 


± 


120 


Apr. 


2.6 ±0.8 


29.1 ±0.7 


11804 ±63 


5028 


± 


155 


May. 


4.6 ± 0.7 


21.2 ±0.6 


11729 ±85 


4465 


± 


181 


Jun. 


6.5 ± 0.4 


14.7 ±0.5 


12107 ± 109 


3358 


± 


167 


Jul. 


5.4 ±0.2 


6.4 ± 0.3 


11852 ± 128 


2933 


± 


184 


Aug. 


5.0 ±0.2 


7.8 ±0.3 


11955 ± 101 


2779 


± 


141 


Sep. 


5.7 ±0.3 


16.3 ±0.4 


12076 ± 80 


3282 


± 


120 


Oct. 


4.8 ±0.5 


22.4 ± 0.5 


11889 ±76 


4158 


± 


144 


Nov. 


0.5 ±0.8 


30.1 ±0.8 


11453 ±50 


5297 


± 


152 


Dec. 


0.0 ± 0.7 


32.3 ±0.7 


11462 ±41 


5295 


± 


123 



Table 6. Flagstaff Wind Speed Model Coefficients 



Month 


AO 


Al 


A2 




A3 


Jan. 


-8.1 ± 4.0 


36.8 ± 3.9 


11486 


± 


63 


6939 ± 528 


Feb. 


1.3 ± 1.6 


29.2 ± 1.5 


11404 


± 


55 


5210 ± 244 


Mar. 


-0.6 ± 0.9 


27.0 ± 0.9 


11487 


± 


29 


5702 ± 155 


Apr. 


2.9 ±0.8 


24.6 ±0.8 


11280 


± 


34 


5106 ± 151 


May. 


6.3 ± 0.7 


17.7 ± 0.7 


11569 


± 


68 


4266 ± 192 


Jun. 


8.7 ±0.3 


11.8 ±0.4 


12084 


± 


80 


2865 ± 123 


Jul. 


6.0 ± 0.2 


7.7 ±0.2 


11847 


± 


68 


2606 ± 101 


Aug. 


5.7 ±0.2 


10.3 ± 0.3 


11967 


± 


61 


2583 ± 88 


Sep. 


8.0 ± 0.4 


16.2 ±0.5 


12097 


± 


82 


3133 ± 137 


Oct. 


4.5 ±0.8 


19.8 ± 0.7 


11669 


± 


59 


4587 ± 183 


Nov. 


-0.2 ± 1.7 


29.0 ± 1.6 


11443 


± 


54 


5484 ± 266 


Dec. 


-2.9 ± 2.8 


31.7 ± 2.7 


11246 


± 


61 


5969 ± 398 



Table 7. Tenerife Wind Speed Model Coefficients 



Month 


AO 


Al 


A2 


A3 




Jan. 


2.5 ± 0.9 


20.3 ± 0.9 


12371 ± 159 


7964 


± 


341 


Feb. 


4.0 ± 0.6 


29.7 ± 0.6 


12102 ± 82 


6478 


± 


149 


Mar. 


1.7 ± 1.1 


27.5 ± 1.0 


12847 ± 173 


8041 


± 


325 


Apr. 


3.0 ± 0.9 


30.6 ±0.8 


12332 ± 120 


6784 


± 


216 


May. 


4.6 ±0.5 


24.9 ± 0.5 


12294 ± 92 


5646 


± 


143 


Jun. 


4.9 ±0.8 


14.9 ±0.8 


12711 ±282 


5899 


± 


420 


Jul. 


6.7 ± 0.4 


8.1 ±0.6 


12536 ± 232 


3760 


± 


304 


Aug. 


6.3 ± 0.4 


8.3 ± 0.7 


12137 ±226 


3630 


± 


296 


Sep. 


6.0 ± 0.4 


11.5 ±0.6 


12137 ± 172 


4341 


± 


241 


Oct. 


-0.2 ± 2.4 


22.9 ± 2.2 


11545 ±218 


8465 


± 


661 


Nov. 


4.5 ± 0.7 


16.2 ±0.7 


12313 ± 191 


6579 


± 


336 


Dec. 


5.6 ± 0.7 


20.9 ± 0.7 


12074 ± 139 


6033 


± 


238 



Table 8. Antofagasta Wind Speed Model Coefficients 



Month AO Al A2 A3 



Jan. 
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13, 
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± 





,3 


12169 


± 


92 


3357 ± 134 


Feb. 
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.5 
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.3 


14. 


,1 


± 





,3 


12090 


± 


86 


3598 ± 131 


Mar. 


4. 


.0 
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directions often show sharp spikes at the edges of jets. Those are an artificial 
result of the averaging of directions and can be ignored. 

For both hemispheres, the tropospheric winds increase in the winter months 
and decrease in the summer months. Stratospheric winds increase in the sum- 
mer months, the opposite of the tropospheric winds. During the summer 
months, Tenerife has a low level (~ 5 km) jet form. This is most obvious in 
July and August, but first appears in June and persists to October. In the Green- 
wood model, Aq is the low altitude wind speed and that is true for our model for 
Hawaii. It is not always true for other sites, as these values are sometimes neg- 
ative. The reason for this is that the wind profiles for those sites has a slightly 
different shape than the Hawaiian wind profile. This can be seen in Figs. [Handle 

Greenwood [9] only provides a model for the wind speed, though ifTOl does 
show the annual mean wind direction for Hilo and Lihue. As shown by Figs. 
[2H8]the wind direction is not an easily modeled function. It can most easily be 
modeled as a square wave. Since stability of flows does depend upon the amount 
of shear present, an analysis that computed the wind shear at the upper and 
lower boundaries would tend to overestimate the instability. Looking at Figs. 
[2HHI the wind direction does not change instantaneously, but it does change 
very quickly, often changing direction by 180° in a km. Some regions of slowly 
changing direction, most commonly above 25 km, are not very well modeled 
by a square wave. A square wave model will suffice for modeling where phase 
screens at only a few altitudes are used to model the distribution of turbulence. 

Note that the reversal of wind direction at altitudes above 20 km is a well 
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Fig. 2. The wind speed and direction for Hilo, Hawaii. Each of the graphs 
in this figure has the the measured wind speed as a function of height plot- 
ted with a solid black line and the wind speed computed from the model 
coefficients in Table |2]is shown as a black dashed line. The wind direction 
is plotted as a solid red line. The wind speed is marked along the bottom 
of the graph, and the wind direction is marked on the top of the graph. The 
first letter of the corresponding month is in the upper right corner. 
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Fig. 3. The wind speed and direction for Oakland, California. The layout 
of the figure is the same as in Fig. [2 except the fitted wind model comes 
from Table |3] 
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Fig. 4. The wind speed and direction for San Diego, California. The layout 
of the figure is the same as in Fig. [2 except the fitted wind model comes 
from Table |4j 
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Fig. 5. The wind speed and direction for Tucson, Arizona. The layout of 
the figure is the same as in Fig. [2 except the fitted wind model comes from 
Table [5j 
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Fig. 6. The wind speed and direction for Flagstaff, Arizona. The layout of 
the figure is the same as in Fig. [2 except the fitted wind model comes from 
Table gj 
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Fig. 7. The wind speed and direction for Tenerife, Canary Islands, Spain, 
which is useful for the observatories located in the Canary Islands. The 
layout of the figure is the same as in Fig. [2j except the fitted wind model 
comes from Table [7] 
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Fig. 8. The wind speed and direction for Antofagasta, Chile. The layout of 
the figure is the same as in Fig. [2 except the fitted wind model comes from 
Table [8J 



known feature of the zonal wind flow. It is due to a higher altitude jet whose 
center is north of the equator in winter, and south of the equator in summer. 
There are also northern and southern jets at higher altitudes which move in 
the same direction as the jets in the lower stratosphere. These are seasonal and 
impact the wind speed and directions in our data lfl4l . 

3. Richardson Number 

The gradient Richardson number is a dimensionless ratio, Ri, related to the 
buoyant production or consumption of turbulence divided by the shear produc- 
tion of turbulence lfT2ll . It is defined as 

Ri= : (4) 

dz) ^ \ dz 

where g is the acceleration of gravity, is the potential temperature, z is the 
altitude, u and v are the horizontal components of the wind vector. Normally 
a gradient Richardson number with a value less than 0.25 is considered to be 
turbulent and a value greater than 0.25 to be non-turbulent. It is used to indicate 
dynamic stability and the formation of turbulence. 

In order to see at which altitudes turbulence forms, we used data from the 
radiosondes to compute the monthly binned Richardson number for each site in 



Table Q] using Eq.@]at each of the fixed grid points used in §[2l The grid heights 
points are every 400m and some of the finer atmospheric layers may not be 
captured with our analysis. The derivatives were computed using a three-point 
Lagrangian interpolation routine. 

The computed Richardson numbers were then binned into two bins: Rj < 
0.25, and Rj > 0.25. For a given site, at each altitude the numbers in each bin 
were normalized by the total number of data points at that altitude. This gives 
the Richardson number as a function of altitude. These are shown in Figs. [9HT31 
As expected, the Richardson numbers above 20 km are very high and have little 
impact on astronomical seeing; as a result we only plotted the values below 20 
km. It is important to note that these values are for radiosondes launched from 
sites at lower altitudes than astronomical observatories. The highest Richard- 
son number values occur at these low altitudes and are not a concern to most 
astronomical observatories. 
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Fig. 9. The binned Richardson number for Hilo, Hawaii The figure shows 
the fraction of Richardson numbers falling into two bins. Yellow signifies 
turbulent conditions (Rj < 0.25) and blue signifies non-turbulent conditions 
(Rj > 0.25). 







Examining the individual months reveals that the figures are not as smooth 
as the wind speed figures, but instead have many fine layers. The layers are 
significant, even if they do not seem to explicitly show turbulence (ie their value 
of Ri is i 0.25). Because the Rj calculation has a height increment of 400m, and 
because the precise positions and altitudes of jets vary over periods of a day 
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Fig. 10. The binned Richardson number for Oakland, California. Yellow 
signifies turbulent conditions (/?,• < 0.25) and blue signifies non-turbulent 
conditions (7?,- > 0.25). 




Fig. 11. The binned Richardson number for San Diego, California. Yellow 
signifies turbulent conditions (/?,• < 0.25) and blue signifies non-turbulent 
conditions (7?,- > 0.25). 
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Fig. 12. The binned Richardson number for Tucson, Arizona. Yellow sig- 
nifies turbulent conditions (/?,• < 0.25) and blue signifies non-turbulent con- 
ditions (Rj > 0.25). 
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Fig. 13. The binned Richardson number for Flagstaff, Arizona. Yellow 
signifies turbulent conditions (/?,• < 0.25) and blue signifies non-turbulent 
conditions (7?,- > 0.25). 
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Fig. 14. The binned Richardson number for Tenerife, Canary Islands, 
Spain. Yellow signifies turbulent conditions (Rj < 0.25) and blue signifies 
non-turbulent conditions > 0.25). 




Fig. 15. The binned Richardson number for Antofagasta, Chile. Yellow 
signifies turbulent conditions (/?,• < 0.25) and blue signifies non-turbulent 
conditions (7?,- > 0.25). 



or less, the effective values of the velocity derivatives in Equation |4] are lower 
in our results than they might be at any given time. This averaging leads to 
relatively large calculated Rj values. It also means that our data will not show 
many instances of Ri < 0.25. But where is does occur it is significant. Moreover, 
any drop in the plotted Rj number indicates a persistent region of instability, and 
it is likely that some turbulence is present there for extended periods, 

Since low gradient Richardson numbers indicate the presence of atmospheric 
turbulence, they can be an indicator of seeing. In Fig. [9] the winter months for 
Hawaii have the highest proportion of Richardson numbers lower than 0.25 
(yellow values), and the summer months have the lowest. In other words the 
winter months have more turbulent conditions than the summer months. Studies 
of seeing from Mauna Kea show that this seasonal behavior of seeing correlates 
with this very well lfT51[T6l[T7ll . There is relatively little seasonal variation in for 
Antofagasta, Chile (Fig.[15]) and this correlates well with seeing measurements 
carried at nearby Cerro Armazones ifTTl which show no clear evidence for any 
seasonal dependency. The trend in Fig. [14] agrees with seeing measurements for 
La Palma, Canary Islands, showing that seeing is best May-September lfT8l[T9ll . 
Flagstaff, Arizona (Fig. [13] has more turbulent values than Oakland, Califor- 
nia, which is agreement with site comparison done between Anderson Mesa 
(outside of Flagstaff) and Lick Observatory (near Oakland) [20J. This shows 
that Richardson numbers computed from radiosonde data can be used to dif- 
ferentiate between astronomical sites at a macro level which can be of use for 
narrowing site selection for new telescopes. It is important to remember that 
the Richardson numbers computed from radiosonde data do not include local 
effects such as the turbulent surface boundary layer that can have a profound 
impact on seeing. Small differences in altitudes and location can alter the ob- 
served seeing significantly |[20l . 

While the Richardson number does not indicate the strength of turbulence, 
knowing that there are persistent turbulent layers may be useful to site planners. 
For example, knowing roughly where and when one might expect upper layers 
of turbulence to be present could be useful for measurements, such as laser 
communications, which are strongly influenced by scintillation. On the other 
hand, we have seen that some layers closer to the ground have R, ;j 0.25. These 
must also be persistent to show up over the averaging. Knowing that these lay- 
ers exist is useful for large telescopes that plan to compensate for near ground 
effects. 

4. Other Wind Models 

There have been a few other studies of wind speed for different sites, though 
they have been done for only a limited number of sites and with different 
methodology. Sarazin & Tokovinin |[2T| suggested that the wind speed at the 



200 mbar pressure isobar (V200) is a good indication of atmospheric turbulence 
at a given site (V200 is well monitored for meteorology). The advantage of this 
is that there are models of V200 for the entire planet, while there are only a fixed 
number of radiosondes launched each day. Several groups have used this idea to 
characterize and compare several astronomical sites ll22l l23l l24l including sev- 
eral sites which are included in this study. During the year, the altitude at which 
V200 occurs varies, this makes it difficult to use the values in Eq. \T\ It is useful 
for site comparisons, but it supplements rather than replaces the knowledge of 
the vertical wind profile. 

Garcfa-Lorenzo et al. Il24ll studied V200 for Paranal and La Silla in Chile, 
San Pedro Martir in Mexico, Mauna Kea in Hawaii and La Palma in the Canary 
Islands. The last two observatories are also characterized in this paper. We show 
a seasonal variation in A\ in Tables [2] and [7] which has the highest tropospheric 
wind speed in the spring and the lowest tropospheric winds in the summer. 
Garcfa-Lorenzo et al. showed that the V2oo w ind speeds at Mauna Kea and La 
Palma have this same seasonal behavior. 

5. Conclusions 

We have presented upper-level wind models computed from archival ra- 
diosonde data suitable for the major astronomical observatories in the United 
States, Chile and the Canary Islands. We find that the commonly used Green- 
wood model is not suitable for sites other than Hawaii. In addition it produced 
a single wind profile for the entire year, when it is more correct to use differ- 
ent wind models over the course of the year. We find that, as may be expected 
from geophysical fluid dynamics, that sites with similar latitudes share similar 
wind profiles. In addition we have computed Richardson number profiles from 
the same datasets. Those results indicate the presence of turbulence at different 
altitudes with seasonal variations that seem to agree with variations in seeing. 
These models and results will be of use to modelers of atmospheric turbulence 
and instrument developers. 
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